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1. Introduction 

The fact that the T is the only lepton massive enough to decay into hadrons makes it an excel- 
lent tool to study QCD, both perturbative and non-perturbative, using the precise hadronic T-decay 
data provided by ALEPH The determination of the QCD couphng a^iMT;) [||-|||, which be- 
comes the most precise determination of a^iMz) after QCD running, is an excellent example. In 
this particular case the non-perturbative contributions are strongly suppressed, but in other analy- 
ses the non-perturbative effects are sizable and then one can extract important phenomenological 
hadronic matrix elements and other non-perturbative QCD quantities. Thanks to the fact that the 
spectral function of the T decay can be separated experimentally in its vector and axial-vector 
contributions, we can study their difference that is specially interesting because it vanishes in per- 
turbative QCD (in the chiral limit) and therefore it is a purely non-perturbative quantity. 

The T-decay measurement of this V — A spectral function has been used to perform [^-^ 



phenomenological tests of the so-called Weinberg sum rules (WSRs) [10], to compute the elec 



tromagnetic mass difference of the pions [pO, and to determine several QCD vacuum condensates 



[11, 12] relevant for the computation of ej^/g^f [13]. The common idea under these studies is the 



use of the analyticity properties of the different two-point correlation functions appearing in the 
dynamical description of the T hadronic width. As it is well known, analyticity allows us to relate 
different regions of the ^^-complex plane. Roughly speaking, one can relate in this way regions 
where we are able to compute analytically, either with Chiral Perturbation Theory (;t;PT) or with 
the short-distance Operator Product Expansion (OPE), with regions where we are not able to com- 
pute (except perhaps in the lattice) but that are experimentally accessible. This connection can be 
used either to predict observables that we are not able to calculate "directly" or, in the other way 
around, to extract the value of QCD parameters that are not fixed theoretically. 



Using ;^PT 1|14|-|16|], the hadronic T-decay data can also be related to order parameters of the 
spontaneous chiral symmetry breaking (S;^SB) of QCD. ;^PT is the effective field theory of QCD 
at very low energies that describes the physics of the S;^SB Nambu-Goldstone bosons through an 
expansion in external momenta and quark masses, with coefficients that are order parameters of 
S;^SB. At lowest order (LO), i.e. ff{p^), all low-energy observables are described in terms of 
the pion decay constant ~ 92.4 MeV and the light quark condensate. At ^{p^), the SU(3) 
;^PT Lagrangian contains 12 low-energy constants (LECs), L,=i ... lo and H12 [jl^, whereas at 
G{jp^) we have 94 (23) additional parameters C,=i ... 94 (Cj^j ... 23) in the even (odd) intrinsic parity 



sector 018|]. These LECs are not fixed by symmetry requirements alone and have to be determined 



phenomenologically or using non-perturbative techniques. Values for the L, couplings have been 



obtained in the past with an acceptable accuracy (a recent compilation can be found in ref. []190), 
but much less well determined are the ^{p^) couplings C,. 

There has been a lot of recent activity to determine the chiral LECs analytically, using as much 



as possible QCD information [g0|-g9|], and from lattice simulations [^-^. This strong effort 
is motivated by the precision required in present phenomenological applications, which makes 
necessary to include corrections of ^{p^) where the huge number of unknown couplings is the 
major source of theoretical uncertainty. 

Here we explain how the determination of some of these LECs can be improved significantly 
using the most recent experimental data on hadronic T decays [[T]]. In particular we will obtain the 
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most accurate results for the ;^PT couplings Lg, Lio and Csv or equivalently, in the SU (2) ;^PT 



language, Z5, /g and C50 [pOl]. Previous work on Lio using T-decay data can be found in refs. 



11| , pip . Our analysis is the first one which includes the known two-loop ;^PT contributions and, 
therefore, provides also the SU(3) (SU(2)) ff{p^) coupUngs Cg? (C50). 

We will first introduce the sum rule relations that we will use, then we will show our results and 
finally we will compare them with other recent analytic results and hadronic T-data determinations. 



2. Sum rule approach 




Figure 1: Analytic structure of Tl{s). 



The basic objects of the theoretical analysis are the two-point correlation functions of the non- 



strange vector 



ud 



ud,^ 



Kd 



u y^d) and axial- vector {^^^^ = A^^ = u J^ysd) quark currents: 



d4^e'-'?Mo|r(^j^,(x)X.(o)0lo) 



(0) 

ud,f 



(2.1) 



(0+1) 



(0+1) 



In particular, we are interested in the difference Tl{s) = H^'^^'^ — H^^^'^ and we will work in the 

isospin limit (wij, = ma) where n[^'y(g^) = 0. The analytic behaviour of this correlator is shown 
in Fig.|l|, together with the complex circuit that we will use to apply Cauchy's theorem. As we 
are interested in relating the ;^PT domain (very low energies) with the T data, we multiply this 
correlator by a weight function of the form I/5" with n> Q. In this way we generate a residue at 
5 = 0. Taking into account the OPE associated with our correlator at large momenta and working 
with the cases n = 1,2, one gets the following sum rules (see ref. [30| for a careful derivation): 
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(2.2) 
(2.3) 



where the integrations start at the threshold Sth = 4m|. These two relations represent the starting 
point of our work and define the effective parameters Ljq and Cfj. Their interest stems from the 
fact that the l.h.s. can be extracted from the data (see Section ^), while the r.h.s. can be rigourously 
calculated within ;^PT in terms of the LECs that we want to determine (see Section H). 
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3. Determination of the effective parameters 




Figure 2: L^g (io) and C|y (io) from different sum rules. For clarity, we do not include the error bands 
associated with the modified weights. 



We use the recent ALEPH data on hadronic T decays [pt], that provide the most precise mea- 



surement of the V — A spectral function. In the integrals of equations (2.2) and (2.3) we are forced 
to cut the integration at a finite value ^o, neglecting in this way the rest of the integral from sq to 
infinity. The superconvergence properties of Tl{s) at large momenta imply a tiny contribution from 
the neglected range of integration, provided ^'o is large enough. Nevertheless, this generates a the- 
oretical error called quark-hadron duality violation (DV)^. From the so-sensitivity of the effective 
parameters one can assess the size of this error. 

In Fig. ^ we plot the value of L^q obtained for different values of (solid lines), with the 
one-sigma experimental error band, and we can see a quite stable result at i^o ^ 2 GeV^. The weight 
function l/s decreases the impact of the high-energy region, minimising the DV. The resulting 
integral is then much better behaved than the corresponding sum rules with s" (n > 0) weights. 

There are some possible strategies to estimate the value of Ljq and his error. One is to give 
the predictions fixing sq at the so-called "duality points", two points where the first and second 
WSRs [|10|| happen to be satisfied. In this way we get Ljf = -(6.50 ±0.13) • 10"^, where the 



uncertainty covers the values obtained at the two "duality points". If we assume that the integral 



(2.2) oscillates around his asymptotic value with decreasing oscillations and we perform an average 
between the maxima and minima of the oscillations we get Ljq = —(6.5 ±0.2) • lO^-'. Another 
way of estimating the DV uses appropriate oscillating functions defined in [ ^ ] which mimic the 
real quark-hadron oscillations above the data. These functions are defined such that they match 
the data at ~ 3 GeV^, go to zero with decreasing oscillations and satisfy the two WSRs. We 
find in this way Ljq = —(6.50 ±0.12) • 10^^, where the error spans the range generated by the 
different functions used. These estimates are in good agreement with each other and give us a first 
determination, but the most precise way to evaluate the error can be obtained taking advantage 
of the WSRs to construct modified sum rules with weight factors w{s) proportional to (1 — s/sq), 
in order to suppress numerically the role of the suspect region around s sq [||]. Fig. ^ shows 
the results obtained with wi{s) = {1—s/sq) /s (dashed line) and W2{s) = {\—s/sq)^/s (dot-dashed 
line). These weights give rise to very stable results over a quite wide range of sq values. One gets 



'From a different but equivalent perspective, we are assuming that the OPE is a good approximation for n(i) at any 
|i|=.?Oi what is not expected to happen near the real axis, and that produces the DV. 
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L^ff = -(6.51 ±0.06) • 10"^ using wi{s) and = -(6.45 ±0.06) • lO^^ using W2is). Taking into 
account all the previous discussion, we quote as our final result: 



L^ff =- (6.48 ± 0.06) -10"^ (3.1) 

We have made a completely analogous analysis to determine C^j. The results are shown in 
Fig. ^. The solid lines, obtained from Eq. (23), are much more stable than the corresponding results 



for LjQ, due to the l/s factor in the integrand. The dashed and dot-dashed lines have been obtained 



with the modified weights W3(s) =p-(^l — and W4(s) =^ V^~'~^^/' "^^^ agreement 

among the different estimates is quite remarkable, and our final result is 

C^^ = (8.18±0.14)-10"^GeV-2_ (-32) 



Our result for Ljq agrees with [S, ITTl 511], but our estimation includes a more careful assessment of 



the theoretical errors. The 3.2a discrepancy between the estimation of ref. ||9|] and ours is caused 
by an underestimation of the systematic error associated with the duality-point approach used in 



that reference. Only in ref. [ ]31[ ] Cgy is also determined and it is in good agreement with our result. 



4. Determination of the ;tPT couplings 



Using the results of ref. [g^] to calculate in ;^PT the r.h.s. of equations ( |2.2| ) and ( |2.3[ ), we get 

-SLl^ = -SUM + GtL{^)+ Gf,L{^) + GU^) + GU^) , (4.1) 
16Qf = HtL+ \6CM+Htt^{pi)+H^^{pi) , (4.2) 

where the functions G\"^{iJ.),H'"^{iJ.) are corrections of order p'" generated at the «-loop level, which 
explicit analytic form | |30| ] is omitted for simplicity. 

Working at ^(p^), the determination of the chiral coupling Lio is straightforward. One gets 

L[q{h =Mp) = -(5.22 ±0.06) • 10"^ . (4.3) 

At order p^, the numerical relation is more involved because the small corrections Gq^ il(m) 
contain some LECs that represent the main source of uncertainty for L\q. It is useful to clas- 
sify the G{p^) contributions through their ordering within the Xj'Nc expansion. The tree-level 
term Go^(/i) contains the only ff{p^) correction in the \wg&-Nc limit, 4m|(Cgj — C[2— Cgg); 
this correction is numerically small because of the m\ suppression and can be estimated with a 



moderate accuracy [^, |32|, |34|, |35|]. At NLO Go^(/x) contributes with a term of the form 
m^(Cg2— C^g— Cgj). In the absence of information about these LECs we will adopt the conserva- 
tive range |Cg2— ^[3- CgJ < |C^j — Cf2— Cgol/3, which generates the uncertainty that will dominate 
our final error on Ljq. Also at this order in \/Nc there is the one-loop correction G\^{jjL), that is 
proportional to Lg which is better known [06|]. Calculating the suppressed two-loop function 



G2i{n) and taking all these contributions into account we finally get the wanted ff{p^) result: 

L\o{Mp) = -(4.06±0.04^eff±0.39LECs)-10"^ = -(4.06 ± 0.39) • 10"^ (4.4) 
where the error has been split into its two main components. 
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A recent reanalysis of the decay —>■ e^vj [p5[], using new experimental data, has provided 
quite accurate values for the combination L9 + L10, both at order and p^, that combined with our 
results forLjQ(Mp) give us 

f (6.54 ±0.15) -10-3 [^(/)], 
L5(Mp) = <^ (4.5) 
[ (5.50 ±0.40) -10-3 [^(Z' )]• 

Repeating the process we have done for Ljq with Qy (where the only LEC involved is Lg) we get 

q,(Mp) = (4.89 ±0.19) • lO^^GeV^^ (4.6) 



5. SU(2) zPT 



Up to now, we have discussed the LECs of the usual SU(3) ;^PT (X^i)- It turns useful to con- 
sider also the effective low-energy theory with only two flavours of light quarks (xPT^i)- In some 
cases, this allows to perform high-accuracy phenomenological determinations of the corresponding 
LECs at NLO. Moreover, recent lattice simulations ^] with two dynamical quarks are already 
able to obtain the SU(2) LECs with sufficient accuracy and this is an important check for them. 

In SU(2) xPT, there are ten LECs, /,=i,..7 and hi^2,3, at ff{p'^) (NLO) |jT5|l. Using the 
relation between l^{p.) and Ljq(/i), recently obtained in ref. [|^], we get 



13.30±0.11 
12.24 ±0.21 



(5.1) 



Analogously to the SU (3) case, the combination Zg — I5 has been determined from the analysis of 
71 — > / V7 [44]. In combination with our determinations for this gives us^ 



15.80 ±0.29 
15.22 ±0.39 



(5.2) 



Making use of the recent results obtained in reference [ p7| ] we can also rewrite our result for 
Cgy in the X^T^2 language, getting in this way the first determination of Cjq 



c^o(Mp) = (4.95 ±0.19) • 10"^ GeV" 



(5.3) 



6. Summary and comparison with previous estimates 

Tables |l| and ^ summarize our determinations of chiral LECs at ^{p^) and ^{p^), respectively. 
They have been obtained through a sum rule analysis that only uses general properties of QCD and 
the measured V —A spectral function [jl]], and taking into account the results of refs. [ p5| , |44] ]. 

Our determination of Ljq is the first one extracted from T-decay data at &{p^). We can 
make an indirect and interesting check comparing our G{p^) result for Lg (Zg) with the value 

^Actually, at order p^, the most precise value of the combination 7(, —J^ is obtained if we calculate it from the SU (3) 
combination L9 of ref. [|33]. In this way we have obtained a prediction for Ig that supersedes that of ref. [BOp. 
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ZPT2 


ZPT3 


I5 = 12.24±0.21 


L'jq(Mp) = -(4.06 ±0.39) • 10-3 


76 = 15.22 ±0.39 


L5(Mp) = (5.50±0.40)-10-3 


4o = (4.95 ±0.19) • 10-3 


C'^j{Mp) = (4.89 ±0.19) • 10-3 GeV-2 



Table 1: Results for the ;^PT LECs obtained at ff{p^). 



ZPT2 


ZPT3 


I5 = 13.30±0.11 


L\q{Mp) = -(5.22 ±0.06) • 10-3 


76 = 15.80 ±0.29 


L^(Mp) = (6.54±0.15)-10-3 



Table 2: Results for the ^PT LECs obtained at ^(p"*). 



L5(Mp) = (5.93 ± 0.43) • 10-3 = 16.0 ± 0.5 ± 0.7) obtained from the charge radius of the 
pion [ |36| ] ([p5|]). The agreement is very good and the improvement in the numerical value of 76 is 
remarkable. 

At order p'^ we do have a previous estimate of L'^q from T data that found L\q{Mp) = 
— (5.13 ±0.19) • 10-3, through a simultaneous fit of this parameter and the OPE corrections of 
dimensions six and eight to several spectral moments of the hadronic distribution. This determi- 
nation is in good agreement with our ^{p^) result. Our quoted uncertainty has an smaller ex- 
perimental contribution and includes a better assessment of the theoretical uncertainties. We can 
also perform an indirect check through the comparison of our ff{p^) result for Lg with the value 
Lg{Mp) = (6.9 ± 0.7) • 10-3 obtained from the charge radius of the pion [ 19]. We see again a very 
good agreement and a clear improvement in the precision. 

If we shift now from phenomenology to theory, we can compare our results with those obtained 
from analytical approaches and lattice simulations. Our determinations of L[q{Mp) and Cgj{Mp) 
agree within errors with the large-A'^c estimates based on lowest-meson dominance [^, ^ 
Lio ^ -3/|/(8M^) ^ -5.4 • 10-3 and Cg? ^ lfl/{32M^) ^ 5.3 • 10-3 GeV-^, and with the re- 



sult of ref. [ |28[ ] for Cg?, based on Pade approximants. These predictions however are unable 
to fix the scale dependence which is of higher-order in l/Nc- More recently the resonance chi- 
ral theory Lagrangian [24, ^] has been used to analyse the correlator H{s) at NLO order in the 
l/Nc expansion. Matching the effective field theory description with the short-distance QCD be- 
haviour, both LECs are determined, keeping full control of the /i -dependence. The predicted values 
L^o(Mp) = -(4.4±0.9) • 10-3 and C^j{Mp) 
ment with our results, although less precise. 



(3.6 ± 1.3) • 10-3 GeV-2 are in perfect agree- 
The most recent lattice calculations find the following results (order p^): 



-(5.2±0.5)-10-3 
-(5.7± 1.1 ±0.7) -10-3 



El 



14.9 ±1.2 ±0.7 
11.9±0.7±1.0 



[42] 



[43] 



(6.1) 



They are in good agreement with our determinations (although still far from the phenomenological 



precision), but for the last one that is slightly smaller. As discussed in ref. [43], this is partly due to 
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the deviation of the lattice determination of the pion decay constant from the one. 

Therefore we can conclude that the different analytical approaches and the various lattice 
calculations agree very well with our precise phenomenological values. 
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